Critical temperature for a-particle condensation in asymmetric nuclear matter. 
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The critical temperature for a-particle condensation in nuclear matter with Fermi surface imbal- 
ance between protons and neutrons is determined. The in-medium four-body Schrodinger equation, 
generalizing the Thouless criterion of the BCS transition, is applied using a Hartree-Fock wave func- 
tion for the quartet projected onto zero total momentum in matter with different chemical potentials 
for protons and neutrons. 

PACS numbers: 21.65.Cd, 67.85.Lm 
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I. INTRODUCTION 

Clustering in nuclei and nuclear matter has become an 
increasingly studied subject recently [H4l4j . Since the a- 
particle is the smallest doubly magic nucleus, it is thus 
very stable and excited heavier nuclei may exist in weakly 
bound states consisting of nucleons, pairs, a-particles, 
etc. d 0, El. Recently, excited states of 12 C and 16 0, 
like the Hoyle state in 12 C, have been described as so- 
called "a condensation" states, which can be pictured by 
three or four a-particles occupying the lowest s-state in 
an effective quartet mean-field potential [T], . 

On the other hand, a-clustering in nuclear matter 
has also been investigated as a possible phase in com- 
pact stars or as a precursor of a-clustering in finite 
nuclei. The critical temperature for the a condensa- 
tion has previously been calculated with the momen- 
tum projected Hartree-Fock approximation, yielding re- 
sults consistent with the one by the Faddeev-Yakubovsky 
method [3]. The critical temperature is an important 
quantity with respect to stellar nucleosynthesis and star 
formation [TBI , [la ]. However, in comparison with our 
previous studies, bulk nuclear matter, like the one in 
neutron stars and other compact objects, is unbalanced 
with respect to the particle numbers of the protons and 
neutrons, implying different chemical potentials. In the 
present paper, we investigate the critical temperature in 
asymmetric nuclear matter. 

The applied method is the same as in our previous work 
of Ref. |14j , where the in- medium four-body Schrodinger 
equation is computed employing Hartree-Fock wave func- 
tions under the constraint of total momentum equal to 
zero. While we used four equal single particle wave func- 
tions in the previous paper for symmetric nuclear matter, 
we shall employ different wave functions for protons and 
neutrons in the present case because of the population 
imbalance with the different chemical potentials for pro- 
tons and neutrons. 



In the next section, we shortly will describe the formu- 
lation which we use to compute the critical temperature. 
In Section Hill we show the numerical results for the crit- 
ical temperature as a function of density and chemical 
potential for various population ratios of protons versus 
neutrons. Besides, we will display and discuss the quar- 
tet wave function at the critical temperature. Finally in 
Section [IV1 we summarize. 



II. FORMULATION 

For the asymmetric nuclear matter case, we will closely 
follow the formulation given in our previous paper [l4| for 
symmetric nuclear matter. 

The Hamiltonian is represented by 



H = H + V ^^e lC \ Cl 
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where indices represent momentum, spin, and isospin, 
«i2,i'2' is the antisymmetric two-body interaction ma- 

k 2 

trix, and £{ = o^r- We formally introduce an effective 
mass m* to account approximately for the exchange term 
of the mean field, whereas the direct term is incorporated 
into the chemical potential. However, since in the follow- 
ing we consider mostly only systems at very low density 
we will disregard mean field effects and take for the mass 
the bare one. 

The four-body Schrodinger equation in the medium 
with the eigenvalue E is given by (see 0, [TM1 for 
details) 

£1234^1234 + ^ ^L234;l'2'3'4''01'2'3'4' = -E^1234 (2) 
1'2'3'4' 

where £1234 = Ylt=i £ » an d ^L234 ; i'2'3'4' is of the form 



Vl234;l'2'3'4' — (1 — /l — /2) 2 W 12,1'2"533'^44' 
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/l — fs) 2 W 13,l'3'^22'<544' 



permutation, (3) 



where is the Fermi-Dirac distribution with different 
chemical potentials for protons and neutrons. Otherwise, 
formally the in-medium four-body equation is exactly the 
same as in the symmetric case. The condition of the tran- 
sition to condensation, known as the Thouless criterion, 
is satisfied with E = 2fi p + 2[i n at T — > T c [20] , where \i p 
(fi n ) is the proton (neutron) chemical potential. 

For simplicity, we consider also for proton-neutron im- 
balanced nuclear matter a spin-isospin-independent two 
body interaction. We do not think that this is much less 
justified than in the symmetric case [Tij where this was 
very successful, indeed. We again employ for the four 
body wave function "01234 the Hartree-Fock ansatz, pro- 
jected on zero total momentum: 



"01234 <Pp{kl)(p p (k 2 )tp n (k 3 )ip n (k 4 )xo 

x (27r) 3 c5(fci + k 2 + k 3 + k 4 ) 



(4) 



where ip T (ki) = ip T (\ki\) is the s-wave single particle wave 
functions for protons (t = p) and neutrons (r = n) 1 
respectively. \o is the spin-isospin singlet wave function. 

Substituting the ansatz of Eq. (j4]) into Eq. ((2J, and 
integrating over superfluous variables, we obtain 



where 
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x(ip p (k 2 )) 2 (ip n {k 3 )) 2 (ip n (k 4 )) 2 
x(27r) 3 (5(fc + k 2 + k 3 + k 4 ), 
d 3 k 2 d 3 k 3 d 3 k 4 d 3 k[ d 3 k' 2 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 
x(l-/ P (fc)-/p(fc 2 ))^ 2 ^^ 

x(27r) 3 (5(fc + fc2 -%[ -k' 2 ) 
xip p (k' 1 )ip p (k 2 )ip p (k 2 )((p n (k 3 )) 2 (ip n (k 4 )) 2 
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C p {k) 



(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 
x(l- f P {k) - fn(k3))v ma ,j, 3 

x(2ir) 3 6(k + k 3 -k[ -k' 3 ) 
Xl fip(k' 1 ){ipp(k 2 )) 2 ip n (k 3 )ip n (k 3 )((p n (k 4 )) 2 
x(2Tr) 3 S(k + k 2 + k 3 + k 4 ), (7) 
d 3 k 2 d 3 k 3 d 3 k 4 d 3 k' 3 d 3 k' 4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 

X(l-/n(*3)-/n(*4))%£ il ^ 



x(2Tr) 3 5(k 3 + k 4 ~ k' 3 - k'^) 

x(27r) 3 <5(fc + k 2 + k 3 + k 4 ) 
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d 3 k 2 d 3 k 3 d 3 k 4 d 3 k' 2 d 3 k' 3 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 

X(l-/p(fc 2 )-/n(fc 3 ))^ 2 fc 3 ^ 

x(2n) 3 5{k 2 + k 3 -k' 2 -k' 3 ) 
X(p p (k 2 )ip p (k' 2 )(p n (k 3 )ip n {k 3 )((p n (k 4 )) 2 
x(2TT) 3 5(k + k 2 + k 3 + k 4 ) (8) 

with the symmetric two-body vertex f£ 2 £ 3 p p ■ The cor- 
responding expressions for the neutrons Am B n , and C n 
are obtained by exchanging indices p -f-> n for A p , B p , 
and C p in Eqs. @, and ©. The Fermi distribution 
functions are 



fr(k) 



1 



(t =p,n). 



(9) 



e (&-Mr)/T +1 

Here, comparing with symmetric nuclear matter [lij ]. 
in the imbalanced nuclear matter case, two coupled equa- 
tions are obtained for protons and neutrons. 



III. NUMERICAL CALCULATION 

As seen from Eqs. ([6]) -([8]), since the wave functions 
<Pp,n(k) are mixed up in A Pi n, Bp,n, and C PiTl , Eq. ([5]) 
constitutes in fact two coupled non linear equations to 
be solved self-consistently by iteration. The critical tem- 
perature is derived from the condition 



d 3 k 



(27T) 



d 3 k 



3 <Pp(k)[(Ap(k) + 3C p {k))<p p (k) + 3B p {k)} 



, .„ip n (k)[(A n (k) + 3C n (k))ip n (k) + W n (k) 



= 0. 



(10) 

)] 
(11) 



Given a set of chemical potentials fi p and fi n one can 
solve the two coupled equations (fTO]) and (fTTj) in adjust- 
ing the temperature, to be identified with the critical 
temperature, T c , as the single parameter. 

For the vertex £, £, in Eqs. ([7]) and ©, we again 

take the separable potential of 14j : 



fel k-2 ,k'^ k' 2 



(1*1-1*2) 

Ae 4 



x (2n) 3 S(k 1 +k 2 -k[ -k' 2 ). (12) 



The parameters A and b are adjusted to the binding en- 
ergy (— 28.3MeV) and to the rms radius (1.71fm) of the 
isolated a-particle; A = -992MeVfm 3 and & = 1.43fm _1 . 
As already mentioned, in the symmetric case our proce- 
dure to solve the in-medium four-body equation with our 
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FIG. 1: Critical temperature as functions of the total chemical 
potential /Utotal = + M« (top) and the total free density 
"total (bottom). The density ratio 5 is in Eq. l|13jl. 

ansatz ^ and the separable interaction (IT^j) gave excel- 
lent agreement with a full Faddeev-Yakubovsky solution 
using the Malfhet-Tjon interaction (MT I-III) [2]. We, 
therefore, think that our procedure is valid in the present 
case as well. 

Fig.[T]shows the critical temperature of a condensation 
as a function of the total chemical potential /itotai = + 
fx n . We see that T c decreases as the asymmetry, given by 
the parameter 

g= n "- n r , (13) 

Tin T" Tip 

increases. This is in analogy with the deuteron case (also 
shown) which already had been treated in [2l|, [22j . On 
the other hand it also is interesting to show T c as a func- 
tion of the free density which is 

"total = n P + n n (14) 
f d 3 k 

r ^l, 

nn = 2 J Wf fn{k) > (16) 

where the factor two in front of the integral comes from 
the spin degeneracy, and f p , n (k) is as in Eq. ©. It should 



be emphasized, however, that in the above relation be- 
tween density and chemical potential, the free gas rela- 
tion is used and correlations in the density have been 
neglected. In this sense the dependence of T c on density 
only is indicative, more valid at the higher density side. 
The very low density part where the correlations play a 
more important role shall be treated in a future publica- 
tion. It should, however, be stressed that the dependence 
of T c on the chemical potential as in the upper panel if 
Fig. [TJ stays unaltered. 

The fact that for more asymmetric matter the tran- 
sition temperature decreases, is natural, since as the 
Fermi levels become more and more unequal, the proton- 
neutron correlations will be suppressed. For small <5's, 
i.e., close to the symmetric case, a condensation (quar- 
tering) breaks down at smaller density (smaller chemical 
potential) than deuteron condensation (pairing). This 
effect has already been discussed in our previous work 
for symmetric nuclear matter [Til . [l7j . For large <Ts, 
i.e. strong asymmetries, the behavior is opposite, i.e., 
deuteron condensation breaks down at smaller densities 
than a condensation, because the small binding energy 
of the deuteron can not compensate the difference of the 
chemical potentials. 

More precisely, for small S's, the deuteron with zero 
center of mass momentum is only weakly influenced by 
the density or the total chemical potential as can seen in 
Fig. [TJ However, as S increases, the different chemical po- 
tentials for protons and neutrons very much hinders the 
formation of proton- neutron Cooper pairs in the isoscalar 
channel for rather obvious reasons. The point to make 
here is that because of the much stronger binding per par- 
ticle of the a-particle, the latter is much less influenced 
by the increasing difference of the chemical potentials. 
For the strong asymmetry S = 0.9 in Fig. Q] then finally 
a-particle condensation can exist up to ntotal = 0.02fm~ 3 
(AHotai = 9.3MeV), while the deuteron condensation ex- 
ists only up to ntotai = 0.005fm" 3 (/itotai = 6.0MeV). 

Overall, the behavior of T c is more or less as can be 
expected. We should, however, remark that the criti- 
cal temperature for a-particle condensation stays quite 
high, even for the strongest asymmetry considered here, 
namely S = 0.9. This may be of importance for the pos- 
sibility of a-particle condensation in neutron stars and 
supernovae explosions [23|, [24[ . 

We also show the single particle wave functions of 
protons and neutrons, entering the quartet wave func- 
tion Q , for various ratios of Fermi surface imbalance and 
chemical potentials in Fig. [21 In most cases of Fig. [21 the 
momentum-space wave functions with negative chemical 
potentials are monotonically decreasing whereas the ones 
with positive chemical potentials have a dip at k = 0. 
However, the momentum-space wave functions also de- 
velop a dip at k — even at a negative chemical potential 
as the asymmetry takes on stronger values. This can be 
seen, e.g. for S = 0.5, /Utotal = 8.93MeV, and 6 = 0.9, 
itotai = — H.OMeV in Fig. [H Furthermore, the neutron 
wave function in fc-space with large positive chemical po- 
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FIG. 2: The momentum-space single particle wave functions (pp , (p n at the critical temperature T c as function of k for 5 = 0.0, 
0.5 0.9, and the real-space wave functions (p p , <p n as a function of r for 5 — 0.9 derived from the Fourier transform of <p p , n (k) 
with ifi Py „(r) = J d 3 fte lfc r (pp,n(fc)/(27r) 3 . The top, middle and bottom figures are for /itotai — Mp + ~ — HMeV, ~ O.OMeV, 
and ~ 9.0MeV, respectively. The wave functions are normalized by J d 3 k<fip tn (k) /(27r) 3 = 1. 



tential develops a node. This behavior is similar to the 
wave functions in Ref. [l4|. As shown in Fig. [51 the dis- 
symmetry of proton and neutron wave functions increases 
as 5 increases. As a consequence, the critical temperature 
decreases, and the a condensation breaks down at a more 
dilute density, see Fig. [1] We also present in most right 
figures of Fig. [2] the proton and neutron wave functions 
in real space. In spite of the sometimes strong dissym- 
metry in momentum space, the proton and neutron wave 
functions are relatively more similar to one another in r- 
space. The neutron wave function develops a node as the 
total chemical potential //total = + Mn increases, but 
the negative values of the wave function remain rather 
moderate. 



IV. SUMMARY 



We reported on the critical temperature for the en- 
particle condensation as a function of the density and 
chemical potential in asymmetric nuclear matter. The 
four-body wave function in the medium is calculated with 
Hartree-Fock wave functions projected onto zero total 
momentum, a procedure which was already very success- 
ful in the symmetric case. Not unexpectedly the tran- 
sition temperature decreases with increasing asymmetry. 
However, it was shown that T c stays relatively high for 
very strong asymmetries, a fact of importance in the as- 
trophysical context. The single particle wave functions 
of proton and neutron were also shown and discussed. 
The neutron wave function in momentum space develops 
a node for strong asymmetries and high densities, a fact 
familiar from ordinary pairing. It also was shown that 
asymmetry affects deuteron pairing more strongly than 
a-particle condensation. Therefore, at high asymmetries, 
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if at all, a-particle condensation dominates over pairing 
at all possible densities. 

The condition of Eqs. (fT0|) and corresponds to 
calculating the critical temperature at the point where 
the a-particles in nuclear matter dissociate into four free 
particles, i.e. at the Mott transition temperature. In the 
higher densities considered here, this temperature is con- 
sistent with the phase transition temperature between su- 
perfluid and normal phase. In the strong coupling limit, 
this is not the case, since the Bose Einstein condensation 
breaks down at lower temperature than the one of disso- 
ciation of a bound state to free particles. This problem 



in the BCS-BEC crossover can be solved by the Nozieres 
and Schmitt-Rink theory, where the pair fluctuations are 
systematically included beyond the mean field approxi- 
mation [25| . The extension of this theory to the quartet 
condensation is difficult but will be attempted in future 
work. 
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